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ABSTRACT 

We present a technique to constrain galaxy cluster pressure profiles by jointly fitting Sunyaev- 
Zerdovich effect (SZE) data obtained with MUSTANG and Bolocam for the clusters Abell 1835 and 
MACS0647. Bolocam and MUSTANG probe different angular scales and are thus highly comple¬ 
mentary. We find that the addition of the high resolution MUSTANG data can improve constraints 
on pressure profile parameters relative to those derived solely from Bolocam. In Abell 1835 and 
MACS0647, we find gNFW inner slopes of 7 = 0.361q! 21 an d 7 = 0-381 q.25j respectively when a and 
/3 are constrained to 0.86 and 4.67 respectively. The fitted SZE pressure profiles are in good agreement 
with X-ray derived pressure profiles. 

Subject headings: galaxy clusters: general — galaxy clusters: individual (Abell 1835, MACS 
J0647.7+7015) 


1. INTRODUCTION 

Galaxy clusters are the largest gravitationally bound 
objects in the universe and thus serve as ideal cosmo¬ 
logical probes and astrophysical laboratories. Because 
the formation of galaxy clusters stems from overdensi¬ 
ties of matter and depends on the cosmic composition of 
the universe, one can constrain cosmological parameters 
such as the matter density of the universe, U m , the mat¬ 
ter power spectrum normalization, erg, and the equation 
of state for dark energy density Ua, w (e.g. Carlstrom 
et al. 2002). 

Within a galaxy cluster, the gas in the intracluster 
medium (ICM) constitutes 90% of the baryonic mass 
(Vikhlinin et al. 2006) and is directly observable in the 
X-ray due to bremsstrahlung emission. At millime¬ 
ter and sub-millimeter wavelengths, the ICM is observ¬ 
able via the Sunyaev-Zel’dovich effect (SZE) (Sunyaev 
& Zel’dovich 1972): the inverse Compton scattering of 
cosmic microwave background (CMB) photons off of the 
hot ICM electrons. The thermal SZE is observed as an 
intensity decrement relative to the CMB at wavelengths 
longer than ^1.4 mm (frequencies less than ^220 GHz). 
At longer radio wavelengths, if relativistic electrons are 
present, parts of the ICM may emit synchrotron emis¬ 
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sion. 

In the core of a galaxy cluster, baryonic physics are 
non-negligible and non-trivial. Some observed physical 
processes in the core include shocks and cold fronts (e.g. 
Markevitch & Vikhlinin 2007), sloshing (e.g. Fabian et al. 
2006), and X-ray cavities (McNamara & Nulsen 2007). 
It is also theorized that helium sedimentation should oc¬ 
cur, most noticeably in low redshift, dynamically-relaxed 
clusters (Abramopoulos et al. 1981; Gilfanov & Syunyaev 
1984) and recently the expected helium enhancement via 
sedimentation has been numerically simulated (Peng & 
Nagai 2009). This would result in an offset between X- 
ray and SZE derived pressure profiles. 

At large radii (R > R 500), 11 equilibration timescales 
are longer, accretion is ongoing, and hydrostatic equilib¬ 
rium (HSE) is a poor approximation. Several numerical 
simulations show that the fractional contribution from 
non-thermal pressure increases with radius (Shaw et al. 
2010; Battaglia et al. 2012; Nelson et al. 2014). For all 
three studies, non thermal pressure fractions between 
15% and 30% are found at ( R ~ R 500 ) for redshifts 
0 < z < 1. However, the intermediate radii, between the 
core and outer regions of the galaxy cluster, offer a region 
where self-similar scalings derived from HSE can be used 
to describe simulations and observations (e.g. Kravtsov 
& Borgani 2012). Moreover, both simulations and ob¬ 
servations find low cluster-to-cluster scatter in pressure 
profiles within this intermediate radial range (e.g. Bor¬ 
gani et al. 2004; Nagai et al. 2007; Arnaud et al. 2010; 
Bonamente et al. 2012; Planck Collaboration et al. 2013; 
Sayers et al. 2013). 

While many telescopes capable of making SZE obser¬ 
vations are already operational or are being built, most 
have angular resolutions (full width, half maximum - 
FWHM) of one areminute or larger. The MUSTANG 
instrument (Dicker et al. 2008) on the 100 meter Robert 
C. Byrd Green Bank Telescope (GBT, Jewell & Prestage 

11 R 500 is the radius at which the enclosed average mass density 
is 500 times the critical density, p c (z ), of the universe 
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2004) with its angular resolution of 9" (FWHM) and sen¬ 
sitivity up to the limit of MUSTANG’S instantaneous 
field of view, 1', is one of only a few SZE instruments 
with sub-arcminute resolution. To probe a wider range 
of scales we complement our MUSTANG data with SZE 
data from Bolocam (Glenn et al. 1998). Bolocam is a 
144-element bolometer array on the Caltech Submillime- 
ter Observatory (CSO) with a beam FWHM of 58" at 
140 GHz and circular FOV with 8' diameter, which is 
well matched to the angular size of R 500 4') for both 

of the clusters in our sample. 

In this paper, we extend the map fitting technique used 
in Young et al. (2014), to simultaneously fit 3D pres¬ 
sure profiles to Bolocam and MUSTANG data. With 
MUSTANG’s high resolution, this is the first analysis to 
make use of SZE observations that cover similar scales 
(O.O 3 R 500 < r < R 500 ) to those probed by X-ray stud¬ 
ies (Nagai et al. 2007; Arnaud et al. 2010, ; hereafter 
N07 and A10 respectively) which have constrained the 
average cluster pressure profile. N07 compared X-ray 
and simulation results over radial scales (O.IR500 < r < 
R500) , whereas A10 used X-ray determined pressure pro¬ 
files for O.O 3 R 500 ;$ r < R 500 ? and simulation results 
for R 500 < r - More recently, the Planck collaboration 
has published an analysis combining XMM-Newton ob¬ 
servations, which span ranges O.O2R500 < r < R 500 with 
Planck observations, which span radial ranges O.IR500 < 
r < 3 R 500 (Planck Collaboration et al. 2013). Addition¬ 
ally, a sample of clusters studied by Bolocam has been 
analyzed using solely the SZE data, which spans radial 
ranges of O.O 7 R 500 < r < 3.bR^oo (Sayers et al. 2013). 

This paper is organized as follows. In Section 2 we 
describe the MUSTANG and Bolocam observations and 
reduction. In Section 3 we detail the method used to 
jointly fit pressure profiles to MUSTANG and Bolocam 
data. We present results from the joint fits in Section 
4 and discuss our results in Section 5. Throughout this 
paper we assume a flat ACDM cosmology with = 
0.3, = 0.7, and H 0 = 70 km s -1 Mpc -1 , consistent 

with the 9-year Wilkinson Microwave Anisotropy Probe 
(WMAP) results reported in Hinshaw et al. (2013). 

2. OBSERVATIONS AND DATA REDUCTION 
2.1. Sample 

To test the application of the joint fitting technique, we 
use two clusters: MACS J0647.4+7015 and Abell 1835. 
Properties of these two clusters are shown in Table 1. 
MACS 0647 (Ebeling et al. 2007, 2010) has been ob¬ 
served by MUSTANG as part of an SZE program to ob¬ 
serve the Cluster Lensing and Supernova with Hubble 
(CLASH) sample (Postman et al. 2012). MACS J0647 is 
well detected by MUSTANG and appears to be relaxed; 
thus it presents a good baseline for testing our fitting 
procedure. We note that it has been previously analyzed 
(Young et al. 2014), albeit with a slightly different ap¬ 
proach. Abell 1835 is a well known cool core cluster (e.g. 
Peterson et al. 2001), with prior SZE detections (e.g. 
Reese et al. 2002; Benson et al. 2004; Bonamente et al. 
2006; Sayers et al. 2011; Mauskopf et al. 2012). MUS¬ 
TANG detects the cluster along with a point source in 
the center of the cluster, thus making it a good cluster to 
demonstrate that the joint fitting technique established 
here can distinguish between point source signal and SZE 


Table 1 

Cluster Properties 


Cluster 2 : 

^500 

-R500 

■M500 


(Mpc) 

0 

(10 14 Mq) 

A1835 0.253 

1.49 ± 0.06 

6.30 

12.3 + 1.4 

MACS J0647.7 0.591 

1.26 + 0.06 

3.16 

10.9 + 1.6 


References. — Mantz et al. (2010) 


signal. 

2.2. MUSTANG Observations 

MUSTANG is a 64 pixel array of Transition Edge Sen¬ 
sor (TES) bolometers arranged in an 8 x 8 array located 
at the Gregorian focus on the 100 m GBT. Operating at 
90 GHz (81-99 GHz), MUSTANG has an angular resolu¬ 
tion of 9" and pixel spacing of 0.63/A resulting in a FOV 
of 42". More detailed information about the instrument 
can be found in Dicker et al. (2008). 



-3-2-10 1 2 3 

Cross Elevation (arcminutes) 


Figure 1. Dotted line: an example GBT trajectory for a 140 
second scan with a Lissajous daisy scan pattern. Three 140 second 
scans result in complete coverage of a circle with 3' radius. The 
plus sign indicates the pointing center used for this daisy scan. 
The diamonds indicate six offset pointings, used on observations 
of MACS 0647. The FOV of MUSTANG is shown by the shaded 
box. The dashed circle encloses 50% of the peak weight over the 
ensemble of scans for a given cluster. 

Observations consisted of scanning on pointing cen¬ 
ters covering the central region of each galaxy cluster. 
A variety of scan durations ranging between 90 and 300 
seconds were employed. Over the seasons of observa¬ 
tions, we found that scans lasting 200 to 220 seconds 
provided the best yields, and accordingly we began to fa¬ 
vor scans of this length. Abell 1835 was observed in the 
winter/spring of 2009 and 2010 with one central pointing 
as described in Korngut et al. (2011). For MACS 0647 
(observed 2011-2013), we adopted an observing strategy 
that consisted of one central pointing, and six off-center 
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Table 2 

Dataset Properties for MUSTANG observations 


Cluster 

Time 

Secondary 

Noise 


(hrs) 

Calibrator 

(/My/beam) 

A1835 

8.6 

1337-1257 

53.3 

MACS J0647.7 

16.4 

0721+7120 

20.8 


pointings, spaced hexagonally such that each was 1' from 
the center, and V away from the other two pointings (Fig¬ 
ure 1). Using the offset pointings with a Lissajous daisy 
scan pattern provides more uniform coverage within the 
central arcminute of our maps than using only one central 
pointing with Lissajous daisy scans. For each pointing 
center, the Lissajous daisy has a 3' radius. The coverage 
(weight) drops to 50% of its peak value at a radius of 
1.3'. 

Absolute flux calibrations are based on the planets 
Mars, Uranus, or Saturn, nebulae, or the star Betelgeuse 
{pt-on)- At least one of these flux calibrators was ob¬ 
served at least once per night. Planets and nebulae are 
the preferred targets, as they can be calibrated directly 
to WMAP observations (Weiland et al. 2011). Betel¬ 
geuse can then be cross calibrated to these planets and 
nebulae, and may be used as an absolute calibration it¬ 
self when no planets or nebulae were observed in a given 
night. We find our calibration is accurate to a 10% RMS 
uncertainty. 

At the start of each night of observing, medium-scale, 
mostly thermal imperfections in the GBT surface are 
measured and corrected using out-of-focus (OOF) holo¬ 
graphic technique (Nikolic et al. 2007). Interspersed with 
scans on clusters, we observe nearby compact quasars 
as secondary calibrators. Secondary calibrators are ob¬ 
served roughly once every 30 minutes, and allow us to 
track the pointing, beam profile, and gain changes of the 
telescope. If the the beam shape or gain degrade by more 
than 10%, another OOF measurement was performed. 
We find that our pointing accuracy is 2". We pair the 
scans on secondary calibrators with off-source scans of an 
internal calibration lamp (CAL) chopped with a 0.5 Hz 
square wave pattern. These CAL scans are taken with 
the telescope at rest. 

All observations of calibrators (primary and secondary) 
also make use of a Lissajous daisy scan pattern. These 
scans have a duration of 90 seconds and radius of 1.5'. 
This scan pattern ensures that each detector passes over 
the point source (within the primary beam) multiple 
times in a scan. Additionally, a radius of 1.5' ensures 
adequate sampling of the error beam. 

Table 2 shows the time spent (on source) observing 
each cluster. MACS 0647, which had exceptionally good 
quality data compared to the other cluster, fares better 
than the expected t -1 / 2 scaling relative to Abell 1835. 
In Table 2 the noise is calculated from a noise map, pro¬ 
duced as described in Section 2.4. Specifically, we smooth 
the map by convolving it with a Gaussian with FWHM 
of 10" and then select the central arcminute (r < 1'), 
where the weight is nearly uniform, and calculate the 
root-mean-square (RMS) of the selected pixels. 

2.3. MUSTANG Beam 


The MUSTANG beam is characterized by making cal¬ 
ibrated maps of compact sources. Point source maps of 
bright secondary calibrators, which are in a high signal- 
to-noise regime, are produced with gentler filtering than 
is used to create galaxy cluster maps(see Section 2.4); the 
primary difference is that no common mode is subtracted 
and we use data from outside the central arcminute to 
calculate and subtract our polynomial fit. For a single 
map, the point source centroid is determined by fitting a 
two dimensional Gaussian to the map using a Levenberg- 
Marquardt least squares minimization in IDL (MPFIT, 
Markwardt 2009). With the centroid determined, pixel 
values and weights are stored as a function of radius, and 
the values are normalized based on the peak of the fit¬ 
ted 2D Gaussian. This constitutes the normalized point 
source profile, or MUSTANG beam profile, for one map. 
2D Gaussians are fit to all 185 secondary calibration 
observations from fall 2011 until spring 2013 with fair 
quality data: wind speeds below 5 m/s, and FWHM 
< 11". The FWHM is taken as the geometric mean of 
the FWHM along the two axes as fit for below. 

To characterize the typical beam for MUSTANG, we 
assume a one dimensional, double Gaussian: 

r 2 r 2 

B(r ) = B\e 2a i + 2ct 2 (1) 

where B\ is the normalization of the primary beam, and 
L >2 is the normalization of the secondary (error) beam. 
This model is then fit, using MPFIT, to the normalized, 
azimuthal profiles of all the secondary calibrators to ob- 

tain B 1 = 0.94+°t *1 = 3.69+1+ B 2 = 0.06±g;^, 

(72 = 12.1±H". This corresponds to a primary beam 
with FWHM of 8.7" and a secondary (error) beam with 
FWHM of 28.4". Figure 2 shows the fitted beam and 
normalized pixel values. 12 The secondary beam is qual¬ 
itatively consistent with the expected near-sidelobes on 
the GBT given the MUSTANG illumination pattern and 
medium-scale aperture phase errors not fully corrected 
by the OOF procedure in Section 2.2. 

Given the dispersion in pixel values in Figure 2, we 
investigate the stability and uncertainty in the fitted pa¬ 
rameters. To investigate the stability of the fitted MUS¬ 
TANG beam, we divide the sample by observing season 
(fall to spring). Fitting a one dimensional, double Gaus¬ 
sian to the beam profile over five seasons of data, we 
find minimal variation and consistent dispersion in pixel 
values. Moreover, the apparent bimodality in the inner 
beam is present in all seasons, and appears to be due 
the ellipticity of individual beams. From the original 
two-dimensional (single Gaussian) fits, we find the mode 
and median of major-to-minor axis ratios are 1.05 and 
1.08 respectively. Taking azimuthal profiles of the fit¬ 
ted (model) two-dimensional Gaussians reveals the same 
bimodality. 

2.4. MUSTANG Reduction 

Processing of MUSTANG data is performed using a 
custom IDL pipeline. Raw data is recorded as time or¬ 
dered data (TOD) from each of the 64 detectors. An 
outline of the data processing for each scan on a galaxy 
cluster is given below. 

12 This is also well approximated as a single Gaussian of FWHM 
9 ". 
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Figure 2. The solid line shows the fitted radial beam profile of 
MUSTANG. Fitting a double Gaussian, the primary beam has a 
FWHM of 8.7", with normalization, B\, of 0.94. The secondary 
beam has a FWHM of 28.4" with a normalization , of 0.06. The 
dashed line shows the weighted RMS of normalized pixels beyond 
27"; it intersects the fitted line at 27". 


(1) We define a pixel mask from the nearest preceding 
CAL scan; unresponsive detectors are masked out. The 
CAL scan provides us with unique gains to be applied to 
each of the responsive detectors. 

(2) A common mode template is calculated as the 
arithmetic mean of the TOD across detectors. The 
pulse tube used to cool the array produces a coher¬ 
ent 1.411 Hz signal across all detectors. A sinusoid 
is used as a template to fit this signal. The common 
mode template, pulse tube template, and a polynomial 
of order N are then simultaneously fit to each detec¬ 
tor and then subtracted. The polynomial order is given 
by N = t scan /t po iy, where t scan is the scan duration, 
tpoly = 2'/f;, and v is the mean scan speed. Typi¬ 
cally, tp 0 iy is roughly 1.0 seconds, while scan durations 
vary between 90 and 180 seconds. Subtracting the com¬ 
mon mode is powerful at removing atmospheric emission, 
but has the downside of removing astronomical signals 
much larger than the instrument FOV. For Abell 1835 
(z = 0.25), 42" corresponds to 166 kpc; for MACS 0647 
(z = 0.59), 42" corresponds to 285 kpc. 

(3) After the common mode and polynomial subtrac¬ 
tion, each scan undergoes further data quality checks: 
spike (glitch) rejection, skewness, and Allan variance. 
Spikes are flagged such that the remaining TOD is still 
used, while detectors that fail skewness and Allan vari¬ 
ance checks are masked for that scan. In total, MACS 
0647 has 11% of its data flagged, and Abell 1835 has 36% 
of its data flagged. 

(4) Individual detector weights are calculated as 1 /of, 
where cq is the RMS of the non-flagged TOD for that 
detector. 

(5) Maps are produced by gridding the TOD in 1" 
pixels in Right Ascension (R.A.) and Declination (Dec). 
A weight map is produced in addition to the signal map. 
Unsmoothed signal-to-noise (SNR) maps are produced 
by dividing the signal map by the inverse square root of 


the weight map. For smoothed SNR maps, the signal 
and variance maps are smoothed, and the SNR map is 
then calculated as the signal map divided by the square 
root of the variance map. 

2.5. MUSTANG Noise 

Because we are in the small signal limit, we need to un¬ 
derstand our noise very well in both the time domain and 
map (spatial) domain. In the map domain, we can pro¬ 
duce noise maps by sending our TOD through the above 
reduction process and either flipping, i.e. reversing, the 
TOD per scan in the time domain, or by flipping sign of 
the gain between scans. The former works by no longer 
allowing the signal to be coherently matched to location 
on the sky, while the latter works by effectively cancelling 
any signal observed. Thorough analysis has shown bet¬ 
ter behavior in the gain-flipped maps. For instance, if 
we make an unsmoothed SNR map from the gain-flipped 
TOD, we find a mean of 0 with a standard deviation of 1. 
In both the gain- and time-flipped noise maps, we find 
that pixel weights accurately reflect the RMS of pixel 
values within a weight bin, and the pixel values follow 
a Gaussian distribution. Smoothed SNR maps produced 
from either of the flipped TOD methods produce stan¬ 
dard deviations greater than 1; the advantage of the gain- 
flipped smoothed SNR maps is that their means are 0, 
whereas the time-flipped SNR maps have means that are 
offset from 0. The standard deviation of our smoothed 
noise SNR maps (as nr) tending towards values greater 
than 1 indicates our smoothing procedure does not ac¬ 
curately handle the weights of each pixel; thus, we use 
&snr to correct our true (not noise) SNR maps. For our 
canonical smoothing kernel (10" FWHM), as nr ~ 1.6; 
we then divide the true SNR map by this factor. As the 
model fitting presented in this work makes use of only the 
non-smoothed maps, this correction factor is only used 
for visualizing smoothed SNR maps. 

2.6. Bolocam Observations and Reduction 

Bolocam is a 144-element camera that was a facility 
instrument on the Caltech Submillimeter Observatory 
(CSO) from 2003 until 2012. Its field of view is 8' in 
diameter, and at 140 GHz it has a resolution of 58" 
FWHM (Glenn et al. (1998); Haig et al. (2004)). The 
clusters were observed with a Lissajous pattern that re¬ 
sults in a tapered coverage dropping to 50% of the peak 
value at a radius of roughly 5', and to 0 at a radius of 10'. 
The Bolocam maps used in this analysis are 14' x 14'. 

The Bolocam data are the same as those used in Cza- 
kon et al. (2014) and Sayers et al. (2013); the details of 
the reduction are given therein, along with Sayers et al. 
(2011). Bolocam observed Abell 1835 for 14.0 hours re¬ 
sulting in a noise of 16.2 [iKcmb -arcminute, and ob¬ 
served MACS 0647 for 11.4 hours resulting in a noise 
of 22.0 IiKcmb~& rcminute. Overall, the reduction and 
calibration is quite similar to that used for MUSTANG, 
and Bolocam achieves a 5% calibration accuracy and 5" 
pointing accuracy. 

3. JOINT MAP FITTING TECHNIQUE 
3.1. Producing Cluster Model Maps 

We choose to construct 3D electron pressure profiles 
as parameterized by a generalized Navarro, Frenk, and 
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White profile (hereafter, gNFW Navarro et al. 1997; Na- 
gai et al. 2007): 


' (C500X)7[l + (C 5 00^) a ]^-^ ^ 

where X = R/R 50 o? and C 500 is the concentration pa¬ 
rameter; one can also write (C 500 X) as (R/R s ), where 
R s = i^soo/Csoo- P is the electron pressure in units of 
the characteristic pressure P 500 • The 3D pressure profile 
is assumed to be spherical and is integrated along the 
line of sight to produce a Compton y profile, given as 

y(r) = ^f T P(r,l)dl (3) 

m e c z J_ OQ 

where R 2 = r 2 + / 2 , r is the projected radius, and l is the 
distance from the center of the cluster along the line of 
sight. Once integrated, y{r) is gridded as y(0) and is re¬ 
alized as a map (pixels of 1" and 20' on a side) using the 
Archive of Chandra Cluster Entropy Profile Tables (AC¬ 
CEPT, Cavagnolo et al. (2009)) centroid for the cluster. 
From here, we produce two model maps: one for Bolo¬ 
cam and one for MUSTANG. In each case, we convolve 
the Compton y map by the appropriate beam shape. For 
Bolocam we use a Gaussian with FWHM = 58", and for 
MUSTANG we use the double Gaussian as found in Sec¬ 
tion 2.3. 

The convolved maps are then regridded to the same 
scale and map size as the reduced data maps for each 
instrument. The regridded Bolocam model map is then 
convolved with a 2D filter function that describes the 
effects of Bolocam data processing (Sayers et al. 2011). 
The MUSTANG map is filtered by converting the model 
map into model TOD, using the true TOD from a galaxy 
cluster as a template (namely for telescope pointing tra¬ 
jectory). The model TOD can then be processed using 
the same custom IDL pipeline used to reduce the data 
to create the filtered MUSTANG model map. 

3.2. Point Source Model maps 

While we do not detect a point source in MACS 0647 
(either in Bolocam or MUSTANG), we clearly detect a 
point source in Abell 1835 in the MUSTANG image. For 
the Bolocam image, the point source in Abell 1835 has 
been subtracted based on an extrapolation of a power 
law fit to the 1.4 GHz NVSS (Condon et al. 1998) and 
30 GHz SZA (Mroczkowski et al. 2009) measurements, 
leading to an assumed flux density of 0.77 =b 0.24 mJy 
(Sayers et al. 2013). 

For MUSTANG, where a point source is detected at 
high significance in our galaxy cluster map, we take the 
following approach: 

1. For each cluster, the same process as in Section 
2.3 is applied to only those secondary calibrators 
which were observed during the same sessions as 
that particular cluster. 

2. The fitted profile is then evaluated as a map, with 
the centroid and total amplitude as determined by 
the 2D Gaussian fit, and shape as determined by 
the ID fit. 

3. The point source map is filtered in the same man¬ 
ner in which the cluster model map is filtered, and 


Table 3 

Parameters of gNFW models 
considered 


Model 

D 500 

a 

p 

7 

Po 

N07 

1.80 

1.30 

4.30 

0.71 

3.94 

A10 

1.18 

1.05 

5.49 

0.31 

7.82 

P12 

1.81 

1.33 

4.13 

0.31 

6.54 

S13 

1.18 

0.86 

3.67 

0.67 

4.29 


Note. — We considered these 
four sets of models and fix a and (3 
for each. 


the resultant image is used as our point source com¬ 
ponent in model fitting, with the normalization as 
a free parameter in the fit. 

3.3. Parameter Space Searched 

Given that the spatial coverage from MUSTANG and 
Bolocam is well suited to constraining the inner pres¬ 
sure profile, we choose to allow the gNFW parameters 7 , 
C 500 , and Pq to vary. To reduce degeneracies, we fix a 
and /3. We choose our fixed parameters from four estab¬ 
lished sets of gNFW parameters: those found in Nagai 
et al. (2007, hereafter N07), Arnaud et al. ( 2010 , here¬ 
after A10), Planck Collaboration et al. (2013, hereafter 
P 12 ), and Sayers et al. (2013, hereafter S13); these are 
summarized in Table 3. We construct a model map for 
each set of <a, /?, 7, and C 500 , and assume a starting 
value for Pq, which is then determined in our fits. Since 
a and /3 are fixed, each cluster is initially searched over 
0 < 7 < 1 in steps of J7 = 0.1, and over 0.1 < C 500 < 2.1 
in steps of JC 500 = 0.1. These ranges are refined after 
the first pass, and generally J 7 is reduced to 0.05. To 
create models in finer steps than J 7 and JC 500 , we inter¬ 
polate filtered model maps from nearest neighbors from 
the grid of original filtered models. 

Given that MUSTANG is sensitive to substructure, we 
fix the MUSTANG cluster model centroid to the position 
of the ACCEPT centroid. If the MUSTANG centroid is 
allowed to vary, we find that the fit can significantly be 
influenced by the cluster substructure. As a result, such 
fits generally do not accurately represent the bulk cluster 
component that we seek to model. However, Bolocam 
maps are dominated by the bulk SZE signal from the 
cluster, and have pointing accuracy of 5". Thus we allow 
the Bolocam pointing is allowed to vary up to a total 
range of 10" in R.A. and Dec relative to the ACCEPT 
centroid with a Gaussian prior of cr cen = 5". 

3.4. Least Squares Fitting 

A set, spanning 7 and C500, of model maps for each of 
the established parameter sets (Section 3.3) is created for 
each galaxy cluster. From each model map, we construct 

a model array of pixel values, ~c^ m , from a subset of the 
pixels in the model maps to fit to the data array of pixel 

values, selected from the same subset of pixels in the 
data map. The subset is chosen to exclude pixels with 
low coverage; principally we select the inner arcminute 
of the MUSTANG maps, and a 14'xl4'box for Bolocam 
data. 

This requires that the data map and model map have 
the exact same astrometry and pixelization. These two 
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arrays can be constructed from the model and data for 
(1) MUSTANG, (2) Bolocam, or (3) MUSTANG and 
Bolocam. In the third case, the arrays are the concate¬ 
nation of the two arrays in the first and second cases. We 
then assume that we can construct our model map as a 
linear combination of components (e.g. a bulk pressure 
profile and in the case of Abell 1835, a point source). To 
do this, we construct an N x M matrix, A, where N is 

the number of data points in our data array, and M 
is the number of model components used. This can be 
written as: 


~d rn, — A CL r. 


(4) 


We can then use the x 2 statistic as our goodness of fit 

X 2 = (7 - 4 m ) r N _1 (7 - itm), (5) 

and find that the minimum y 2 is achieved when: 

= (A r N _1 A) _1 A r N _1 ~<t. (6) 

Here, N is the covariance matrix, which is formally de¬ 
fined as: 


N i:j =< didj > - <di>< dj > (7) 

Here, di is the pixel value in a noise realization map, 
and the average is taken for a given pixel of several noise 
realization maps. 

In practice, we take to be a diagonal matrix for 
MUSTANG and Bolocam, with = Sij/wi , where Sij 
is the Kronecker delta, and Wi is the weight for pixel 
i. MUSTANG’s detector noise is dominated by phonon 
noise; thus with the common mode (which is dominated 
by atmospheric noise) subtracted, we expect the map 
pixel noise to be uncorrelated. Indeed, taking jackknife 
resampling of MUSTANG data has been used to pro¬ 
duce a correlation matrices for clusters and verifies that 
assuming uncorrelated pixel noise in MUSTANG maps 
is valid. 

If we simply append the MUSTANG model array to 
the Bolocam model array: 

m — [3 m, Bolocam m MUSTANG\ (8) 


then solving this set of equations would be a linear prob¬ 
lem. However, we add a calibration offset term, fc, to 
allow for offsets in calibration between MUSTANG and 
Bolocam data, and create the following array: 

m,cal — m, Bolocam •> ^ m, MUSTANG } (9) 


We cast Equation 9 as such because we expect the nor¬ 
malization given by MUSTANG and Bolocam to be 
equal, except for any calibration offset between instru¬ 
ments. Furthermore, we can quantify the calibration un¬ 
certainties and thus put a prior on it. Solving Equation 6 
is no longer a problem of linearly independent variables. 
Thus, we use MPFIT to quickly solve for ~ct m and ob¬ 
tain a pseudo-y 2 (x 2 ). To calculate x 2 , we use the same 
formulation as above, but redefine 


cal — 

m,cal m,old'> A new & m, 



Figure 3. Results from bootstrap Monte Carlo fits to 1000 real¬ 
izations for Abell 1835 with A10 values for a and f3. Here we can 
determine the scatter in each of the fitted gNFW parameters. Sqq 
is the flux density at 90 GHz as determined by MUSTANG. 


where ~ctrn,cai = [~^m, k\, and expand N to allow for the 
extra fitted value. N is modified simply by assigning zero 
for all off-diagonal terms, and adding the expected vari¬ 
ance of k on the diagonal term. We adopt a calibration 
uncertainty of 11.2%, which accounts for the calibration 
accuracies of Bolocam and MUSTANG cited in Section 
2 . 2 . 


3.5. Uncertainties 

For each set of fixed a and (3 (see Table 3), the joint fit 
procedure returns a x 2 . This is used to plot contours of 
Ax 2 , and the equivalent confidence intervals. Following 
Sayers et al. (2011), the uncertainty for each parameter 
is estimated via the lcr confidence interval of the best 
fits over 1000 noise realizations added to model clusters. 
Furthermore, because k is not an independent parameter, 
X 2 will not provide a fully accurate assessment of the 
confidence intervals. 

We create 1000 realizations about the best fit model 
by adding 1000 instances of noise (for both Bolocam and 
MUSTANG) to the best fit model. Then, we find the 
best fit to each realization. The results of these fits is 
shown in Figure 3. The noise realizations for Bolocam 
are precomputed from previous work (Sayers et al. 2013). 
Noise realizations for MUSTANG are computed by ran¬ 
dom number generation with Gaussian distribution for 
each pixel, with cq = 1 /wi , based on the same weight 
used to calculate the noise matrix. These fits are tabu¬ 
lated and used to rescale the Ay 2 confidence intervals as 
seen in Figure 4. The rescaling of the confidence inter¬ 
vals is primarily due to the non-diagonality of the noise 
in the Bolocam covariance matrix. 

We also investigated the potential impact from the un¬ 
certainty in the point source. For a given cluster, we fol¬ 
low the steps in Section 3.2, and we then calculate point 
source uncertainty models which adopt the lcr values for 
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Figure 4. Confidence intervals for Abell 1835 for A10 values of 
a and (3. Top panel shows the confidence intervals for Bolocam 
only fit. Bottom panel shows the confidence intervals for the joint 
fit. Dashed lines show confidence intervals based on Ay 2 . Solid 
lines are the true confidence intervals based on the bootstrap MC 
results. The asterisks denote the best fit. 

the width of the main beam, as reported in Section 2.3, 
and fit the remaining components. The fitting proce¬ 
dure is then rerun twice: once with each of these models. 
Neither the fitted gNFW parameters nor the confidence 
intervals change. However, across the three point source 
models (two uncertainty and best-fit point source mod¬ 
els), the minimum y 2 does change in the expected man¬ 
ner: it is greater for both of the uncertainty models than 
the best fit point source model. 

4. RESULTS 
4.1. Abell 1835 (z=0.25) 

Abell 1835 is a well studied massive cool core clus¬ 
ter. The cool core was noted to have substructure in the 
central 10" by Schmidt et al. (2001), and identified as 
being due the central AGN by McNamara et al. (2006). 
Abell 1835 has also been extensively studied via the SZE 
(Reese et al. 2002; Benson et al. 2004; Bonamente et al. 
2006; Sayers et al. 2011; Mauskopf et al. 2012). The 
models adopted were either beta models or generalized 


beta models, and tend to suggest a shallow slope for the 
pressure interior to 10". Previous analysis of Abell 1835 
with MUSTANG data (Korngut et al. 2011) detected the 
SZE decrement, but not at high significance, which is 
consistent with a featureless, smooth, broad signal. Our 
updated MUSTANG reduction of Abell 1835 is shown 
in Figure 5, and shows the same features as in Korngut 
et al. (2011). 
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Figure 5. MUSTANG SZE specific flux density map of Abell 
1835. The white contours are at the 3 and 4 cr significance level. 
The ACCEPT centroid is the red asterisk. The positive region 
close to the asterisk is the point source. The black circle in the 
upper right represents the effective beam of FWHM of 13.5" (i.e. 
our 9" telescope beam, smoothed by a lO" FWHM Gaussian.) 

We find the best fit A10 model (all parameters but Pq 
fixed to A10 values) has y 2 = 12861. To calculate y 2 
for no cluster model, we fix the point source amplitude 
to that found in the previous fit. We find Ay 2 = 892.1, 
with ADOF= 2 corresponds to 29.7cr significance. For 
Bolocam only, the Ay 2 between a cluster model being fit 
or not yields a 28.9<r detection, while the for MUSTANG 
we find a 10.2cr detection of an A10 model from Ay 2 . 
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Table 4 

Best fit gNFW parameters for Abell 1835 


Table 5 

Abell 1835 point source flux densities 


Model 

P500 

a 

p 

7 

Po 

k 

x 2 

N07 

i-44 tg;S 

1.30 

4.30 

n 74+0.15 
U - '^-0.13 

3-43++ 

1.18 

12837 

A10 

0.83tg;?I 

1.05 

5.49 

0-75++ 

2-54±i;l* 

1.15 

12835 

P12 

1 45+ 0 - 35 

1.33 

4.13 

0.84±g;ii 

2-80±i;fg 

1.14 

12838 

S13 

2 29 +1 ' 30 

0.86 

3.67 

0.36t°;i? 

IQ 0+9-75 
iy, °— 6.16 

1.19 

12831 


Note. — 7 , Po, C 5 oo, and k, the calibration offset, were 
varied. The degrees of freedom were 12880. 


Our best joint fit over the four sets of a and /?, shown 
in Table 4, comes from the S13 values of a and f3\ the 
best fit floated parameters are: 7 = 0.36, Po = 19.3, and 
C 500 = 2.28. Despite variations in the best fit 7 values, 
Figure 6 shows the best joint fit pressure profiles of Abell 
1835 for each of the four sets of fixed a and /?, are in good 
agreement with each other. Moreover, we find that the 
four model fits achieve minimum scatter at two separate 
radii, roughly corresponding to the geometric mean be¬ 
tween the resolution and FOV of each instrument. 

We find the point source in the MUSTANG map at 
R.A.=14:01:02.1, Dec=02:52:47 is best fit with a flux 
density of 1.38 ± 0.10 mJy, and has a correlation coef¬ 
ficient of 0.076 with the cluster amplitude. This minimal 
degeneracy can also be seen in Figure 3. Similarly, chang¬ 
ing the assumed beam shape as discussed in Section 3.5 
has a negligible change on the flux density. The ampli¬ 
tude of the point source suggests a slight flattening of the 
spectral index between v = 1.4 GHz and v = 28.5 GHz 
(Condon et al. (1998),Reese et al. ( 2002 )) of ol v = 0.89, 
to a spectral index of ol v = 0.59 between v = 28.5 GHz 
and v — 90 GHz. Such a spectral index is also consistent 
with McNamara et al. (2014), which find a spectral in¬ 
dex of a u = 0.54 between v = 92 GHz and v = 276 GHz. 
The assumed flux density of the subtracted point source 
for Bolocam, 0.77 =b 0.24 mJy at 140 GHz, is consistent 
with the other measurements. 

The point source flux density found with MUSTANG 
is consistent with that obtained from observations with 
the Atacama Large Millimeter/Sub-millimeter Array 
(ALMA) in McNamara et al. (2014), which find the cen¬ 
tral continuum source has a flux density of 1.26 ±0.03 
mJy. We note they also detect a 1O 1 O M 0 molecu¬ 
lar outflow at 92 GHz, with a total integral flux of 
3.6 Jy km s -1 for CO ( 1 - 0 ), which would correspond 
to an equivalent continuum flux density of 6 fi Jy, and 
would not contribute much additional flux density to the 
point source flux density as seen by McNamara et al. 
(2014). However, this source is reported with a posi¬ 
tion of R.A.=14:01:02.083, Dec=02:52:42.649, which is 
4" offset from the position found in our MUSTANG data. 
Since we consider our positional uncertainty to be 2 ", this 
is a larger than typical pointing offset, but is difficult to 
rule out. A list of selected point source flux densities is 
provided in Table 5. 

Figure 4 shows that Bolocam, by itself, can place fairly 
tight constraints on the gNFW model parameters, pri¬ 
marily on 6500 - That is, for A10 values of a and /?, as in 
Figure 4, Bolocam finds C500 = 0 . 73 to.i 7 , 7 = O^to^, 
and Pq = 1.75l?;g, whereas the joint fit yields C500 = 


Si .4 (mJy) S b 8 5 (mJy) S 90 (mJy) S § 2 (mJy) S § 76 (mJy) 
41.4 ± 1.9 2.76 ±0.14 1.38 ±0.10 1.26 ± 0.03 0.7 ±0.1 

References. — a From Condon et al. (1998) b From Reese 
et al. (2002) c From McNamara et al. (2014). 

Note. — S 90 is from this work. 

0.83to.i55 7 = 0.75^0.17, and Po = 2 . 54 ^ 37 . Across 
the model sets, the trend for Abell 1835 is that joint 
fit tends to loosen the constraint on C 500 , while improv¬ 
ing the constraints on 7 and Po relative to the fits to 
solely Bolocam data. It is worth noting that for each 
value of 7 and C 5005 only Po is a free parameter in the 
Bolocam only fits. In contrast, for the joint fit, the cal¬ 
ibration offset and MUSTANG point source amplitude 
are additional free parameters. The addition of MUS¬ 
TANG data slightly reduces the inner slope, 7, relative 
to the Bolocam-only fit. This is suggestive that Bolo¬ 
cam, with subtraction of its adopted point source model, 
has not underestimated the SZE signal. Moreover, given 
the peak decrement of —20 mJy in the Bolocam map, an 
adjustment of ~ 0.2 mJy, which is the uncertainty on the 
assumed 0.77 mJy, would negligibly alter the constraints. 
Both the Bolocam and joint constraints indicate a rela¬ 
tively steep slope, which is typical for a cool core cluster 
(e.g. Arnaud et al. 2010; Sayers et al. 2013). 


Physical Radius (kpc) 

10 100 1000 



Figure 6. Best fit pressure profiles for Abell 1835 for the different 
sets of fixed a and /3, denoted by the model it is taken from. The 
circles denote radial ranges where the pressure profiles show (local) 
minimum scatter. 


4.2. MACS J0647.7+7015 (z=0.59) 

MACS 0647 is at z = 0.591 and is classified as nei¬ 
ther a cool core nor a disturbed cluster (Sayers et al. 
2013). It was included in the CLASH sample due to its 
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Figure 7. MUSTANG SZE flux density map of MACS 0647. The 
white contours are at the 3, 4, and 5 cr significance level. The 
ACCEPT centroid is the red asterisk. The black circle in the upper 
right represents the effective beam of FWHM of 13.5 X/ . 


Table 6 

Best fit gNFW parameters for MACS 0647 






1. 

1. 

0. 

0 . 

0 . 

0. 

0 . 

1. 

1. 

0 . 

0 . 

0 . 

0 . 

0 . 



Model 

P 500 

a 

0 

7 

Po 

k 

X 2 

N07 

0 93 +0 ' 31 
u.yo_ 0 36 

1.30 

4.30 

0 70+ 0 - 10 
u - ‘ u - 0.17 

2 - 10 ii'.i7 

1.14 

12845 

A10 

0 60 +0 ' 25 
U.OU — 0 22 

1.05 

5.49 

0 61 +0 ' 12 
U - Oi - 0.15 

2 24+ 2,03 
- 1.20 

1.14 

12844 

P12 

1 no + 0.32 
l.OO_0 40 

1.33 

4.13 

0 70+ 0 - 10 
u - ‘ u - 0.17 

9 qc + 1.04 
z ^°_ 1.32 

1.14 

12845 

S13 

1 19+ 0,54 
J " iy -0.64 

0.86 

3.67 

n qo+ 0.20 

U.C)O _0 25 

q i 0 + 4.68 
°- i °- 1.13 

1.14 

12843 


Note. — 7 , Pq , and C 500 were varied. The degrees of freedom 
were 12914. 


strong lensing properties (Postman et al. 2012 ). Gravi¬ 
tational lensing (Zitrin et al. 2011), X-ray surface bright¬ 
ness (Mann & Ebeling 2012 ), and SZE (MUSTANG, see 
Figure 7, and Bolocam) maps all show elongation in an 
east-west direction. In the joint analysis presented here, 
we see that the spherical model provides an adequate fit 
to both datasets and we note that the spherical assump¬ 
tion allows for a easier interpretation of the mass profile 
of the cluster. 

To calculate how significantly Bolocam and MUS¬ 
TANG detect an SZE bulk decrement, we calculate a 
Ay 2 as we did for Abell 1835. Ay 2 is computed as the 
difference of y 2 from an A 10 profile (all parameters but 
Po fixed to A 10 values) fit to our datasets, and x 2 as- 
suming no model. The joint fit (both data sets) yields a 
26.3cr significance, while for Bolocam only, Ay 2 yields a 
23.9cr detection, and in the MUSTANG data we find a 
10 . 8(7 detection of an A 10 model. 

Our best fit model comes from the S13 values of a and 
/?, and finds 7 = 0.38, C 500 = 1.19, and P 0 = 8.18. The 
best joint fits, listed in Table 6 to the four sets of a and /3 


0.5 1.0 1.5 2.0 

C 500 


Figure 8. Confidence intervals for MACS 0647 with A10 values 
of a and (3. Top panel shows the confidence intervals for Bolocam 
only fit. Bottom panel shows the confidence intervals for the joint 
fit. Dashed lines show confidence intervals based on Ay 2 at face 
value. Solid lines are the true confidence intervals based on the 
bootstrap MC results. 


differ by Ay 2 < 3. With Young et al. (2014) constrain- 
ing 7 = 0 . 90 to;° 4 , it might appear that their result is 
significantly discrepant with our best fit 7 = 0 . 611 q!iI 
from the A10 set, even though Young et al. (2014) used 
the identical SZE data as we have used in this analy¬ 
sis. A crucial distinction in the fitting procedures is the 
parameter space searched: in Young et al. (2014), Bolo¬ 
cam is first fit over grid of fixed 7 values, k is fixed at 
1 . 0 , and the parameters C500 and Po are allowed to float. 
The resultant fits for each 7 are then fit as-is (nothing 
allowed to float) to the MUSTANG map. Thus, the re¬ 
ported error bars reflect a one-parameter search, without 
the degeneracies between C500, Po, and 7 folded into it, 
and do not include the y 2 values from the Bolocam fit. 

Given that MUSTANG is only able to constrain the 
pressure profile on scales 9" < 0 < 42", and for MACS 
0647, P 500 = 3.16'and C 500 = 1-18 (A 10 and S13 value) 
or C500 = 1.8 (N07 and P 12 value), then /? should not re¬ 
late to the slope within the scales probed by MUSTANG. 
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It is possible for a to relate to the slope within the scales 
in question. However, As C 500 decrease, especially be¬ 
low 1.0, as is the case in both this work and Young 
et al. (2014), then a will relate less to the slope within 
the scales probed by MUSTANG. Crucially, that cluster 
pressure profile steepen with increasing radius gives rise 
to the degeneracy observed between C 500 and 7 as seen 
in Figures 4 and 8 . 

Figure 8 shows that Bolocam does not place strong 
constraints on 7 and C 500 , especially relative to the joint 
fit (Table 6 ). Specifically, for the A10 set of a and /3, 
the Bolocam finds C 500 = 0.66^o‘38, 7 = 0.39lo‘3g, an d 
Po = 17.7tg;27- A s in Figure 6 , we see in Figure 9 
that the best joint fit pressure profiles from the differ¬ 
ent gNFW fits of MACS 0647 are in good agreement, 
and the radii where Bolocam and MUSTANG have the 
tightest constraints are similar to the radii of tightest 
constraints found in Abell 1835. Figure 10 shows the 
two-dimensional confidence intervals over the parameter 
space searched for MACS 0647. 


Physical Radius (kpc) 
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Figure 9. Best fit pressure profiles for MACS 0647 for the differ¬ 
ent sets of fixed gNFW a and /3, denoted by the model the values 
are taken from. The circles denote radial ranges where the pressure 
profiles show (local) minimum scatter. 


4.3. Comparison with ACCEPT 

ACCEPT (Cavagnolo et al. 2009) utilizes the Chandra 
Data Archive (CDA) to derive entropy (and pressure) 
profiles for 239 galaxy clusters. We summarize how Cav¬ 
agnolo et al. calculate pressure profiles here. ACCEPT 
pressure profiles are calculated as the product of their 
derived temperature and electron density profiles. All 
the ACCEPT data pulled from CDA for ACCEPT were 
taken with the ACIS detectors (Garmire et al. 2003). 
Cavagnolo et al. first derive temperature profiles from 
spectra. The spectra are fitted from annuli containing 
at least 2500 counts after background subtraction. Elec¬ 
tron density is then calculated from surface brightness 



Figure 10. Results from Monte Carlo fits to 1000 realizations for 
MACS 0647 with A10 values for a and f$. 

profiles using the 0.7-2.0 keV range. This time, annuli of 
5" are used (pixels are 0. // 492). To account for the slight 
temperature dependence of X-ray surface brightness, an 
appropriate temperature for each surface brightness an¬ 
nulus is interpolated from the radial temperature profile 
grid. 

We see in Figures 11 and 12 that our joint fits of SZE 
observations agree well with ACCEPT data. We note 
that the ACCEPT profiles show some deviation from a 
gNFW curve. For Abell 1835, there is a bump in the 
pressure at r ^ 40". In MACS 0647, from r ~ 15" to 
r ~ 100", the ACCEPT pressure profile is almost a pure 
power law. 


Physical Radius (kpc) 
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Figure 11. Pressure profiles for Abell 1835. The dotted lines 
show the lcr error bars on the SZE determined pressure profile. 
The best joint fit from this work is the solid line. ACCEPT error 
bars are at the 90% confidence level. 


5. DISCUSSION 
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Figure 12. Pressure profiles for MACS 0647. The best fit (this 
work) is the solid black line, the dashed lines show the la error bars 
from the joint fit, and the solid dark green is the best fit found in 
Young et al. 2014. ACCEPT error bars are at the 90% confidence 
level. 

Here we have presented an approach to jointly fit mod¬ 
els to SZE data from different instruments. With regards 
to our current sample, we have shown that this approach 
finds pressure profiles that are in good agreement with X- 
ray derived pressure profiles. We choose to vary gNFW 
parameters that correspond to the central regions (Po, 
C500, and 7 ) given that we expect our data will be best 
suited to constraining these, and to avoid encountering 
the large degeneracy between all the parameters. 

In Abell 1835 and MACS0647 standard A 10 mod¬ 
els are detected at 29.7 a and 26.3cr significance respec¬ 
tively. For Abell 1835, our best fit is 7 = 0.36^o!36> 

C500 = 2.29to!68> an d Po = 19.3tg° 5 3 with fixed val¬ 
ues of a = 0.86 and /3 = 3.67 from S13. For MACS 
0647, our best fit is 7 = 0.38^0.25> C 500 = 1-19^0.64> anc ^ 
Po = 8 . 18^13 with fixed values of a = 0.86 and (3 = 3.67 
from S13. The error bars have been calculated via 1000 
Monte Carlo realizations and are marginalized over the 
other parameters. While the spread in fitted parameters 
appear large (e.g. in Abell 1835, we find 7 = 0.36 with 
S13 values of a and /3, and we find 7 = 0.84 for P12 val¬ 
ues of a and /?), we note that the fitted pressure profiles 
are in very good agreement (Figures 6 and 9), and the 
profiles reveal nodes of least scatter set by MUSTANG 
and Bolocam. These two nodes occur roughly where ex¬ 
pected: at the geometric mean each instrument’s resolu¬ 
tion and FOV. 

With this approach, we find that contaminants such 
as point sources can be well modeled and their am¬ 
plitudes are minimally degenerate with cluster pressure 
models. The pressure profile constraints do not depend 
on the MUSTANG beam shape assumed to model point 
sources. Furthermore, we find joint fits with MUSTANG 
and Bolocam allow tight constraints on the pressure pro¬ 
files over a large range of scales. 

Confidence intervals for both MACS 0647 and Abell 
1835 show that the constraints are weaker than indicated 


from x 2 . This discrepancy arises because the true covari¬ 
ance matrix for Bolocam is not diagonal, as assumed for 
individual fits. 

Abell 1835 and MACS 0647 also illustrate the impor¬ 
tance of spanning a range of angular scales. While Bolo¬ 
cam data alone provide good constraints on 7 and C 500 
for Abell 1835, Bolocam data do not provide precise con¬ 
straints on these parameters for MACS 0647. While data 
quality and quantity may explain some of the difference, 
a critical distinction is that P 500 occurs at 6.30' for Abell 
1835, and 3.16' for MACS 0647, which is primarily due 
to Abell 1835 being at lower redshift (z = 0.253) than 
MACS 0647 (z = 0.591). The inner pressure profile of 
Abell 1835 can be constrained by Bolocam, but due to 
its larger angular extent, MUSTANG will filter out much 
(r > 1') of the profile within P 500 . However, MACS 0647 
appears smaller on the sky and higher resolution (than 
58") is necessary to constrain the inner pressure profile. 
Therefore, the addition of high resolution MUSTANG 
data becomes crucial for clusters at large redshifts. 

A significant improvement to this approach would be 
to reduce the computational resources required. Cur¬ 
rently, the forward filtering of our MUSTANG models is 
the most computationally expensive step. As this step is 
required for each model created, this restricts how many 
models we can produce, and thus it restricts the param¬ 
eter space we can search. Computing a transfer function 
for MUSTANG would significantly reduce the total com¬ 
putational expense; however, the applicability of such a 
transfer function must be assessed. 

We plan to expand this joint fitting technique to a 
larger sample of clusters, and we note that both Bolocam 
and MUSTANG show good internal consistency in their 
flux calibration. As a result, a single value of k should be 
valid for all of the clusters observed by both instruments, 
and therefore we expect it can be tightly constrained 
using the full sample. This will effectively eliminate one 
of the free parameters in our fits, and produce tighter 
constraints on the gNFW model parameters. 
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